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ABSTRACT 

In  this  paper  the  parameterization  of  an  integrated  electro¬ 
thermal  model  for  an  A123  26650  LiFePO 4  battery  is  presented. 
The  electrical  dynamics  of  the  cell  are  described  by  an  equivalent 
circuit  model.  The  resistances  and  capacitances  of  the  equivalent 
circuit  model  are  identified  at  different  temperatures  and  SOC’s, 
for  charging  and  discharging.  Functions  are  chosen  to  charac¬ 
terize  the  fitted  parameters.  A  two-state  thermal  model  is  used 
to  capture  the  core  and  surface  temperatures  of  the  battery.  The 
electrical  model  is  coupled  with  the  thermal  model  through  heat 
generation.  Parameters  of  the  thermal  model  are  identified  us¬ 
ing  a  least  squares  algorithm.  The  electro-thermal  model  is  then 
validated  against  voltage  and  surface  temperature  measurements 
from  a  realistic  drive  cycle  experiment. 

1  INTRODUCTION 

Lithium  Ion  Batteries  are  attractive  energy  storage  devices 
for  Hybrid  Electric  (HEV),  Plug  In  Hybrid  Electric  (PHEV),  and 
Electric  Vehicles  (EV)  due  to  their  reasonable  power  and  energy 
density.  The  ability  to  accurately  predict  the  electrical  and  tem¬ 
perature  dynamics  of  a  battery  is  critical  for  designing  onboard 
battery  management  systems  (BMS),  and  thermal  management 
systems. 

Electrical  models  vary  in  complexity.  For  some  applications, 
a  simple  model  capturing  the  basic  electrical  behavior  can  be  suf¬ 
ficient  (eg.  an  OCV-R  model).  There  are  more  complex  electro¬ 
chemical  models  [1-3]  that  are  highly  accurate  [4-6],  but  hard 
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to  be  fully  parameterized  [6],  and  require  large  computational 
capacity.  Therefore,  they  are  not  suitable  for  control  oriented 
modeling.  Equivalent  circuit  models  are  commonly  used,  which 
offer  a  tradeoff  between  accuracy  and  simplicity,  and  are  suitable 
for  control  oriented  applications  [7-11]. 

The  equivalent  circuit  model  can  capture  the  terminal  volt¬ 
age  of  the  battery  and  has  been  widely  adopted  since  the  work 
in  [12],  The  voltage  supply  in  the  equivalent  circuit,  shown 
in  Fig.  1,  represents  the  open  circuit  voltage  ( Vocv )  which  is  a 
function  of  state  of  charge.  The  series  resistance  (Rs)  represents 
internal  resistance  of  the  battery.  The  voltage  drop  across  the 
two  resistor-capacitor  (RC)  pairs  (Vi  and  Vf)  are  used  to  model 
the  dynamic  voltage  losses  due  to  lithium  diffusion  in  the  solid 
phase  and  in  the  electrolyte  [13].  These  circuit  elements  depend 
on  state  of  charge  (SOC),  temperature,  and  current  direction  as 
shown  in  [10].  These  parameter  dependencies  are  important  for 
accurately  capturing  the  dynamics  of  battery  terminal  voltage 
throughout  a  usable  range  of  temperature  and  state  of  charge. 

In  addition  to  predicting  the  terminal  voltage,  an  accurate 
model  of  the  battery  temperature  is  needed  for  control  and  ther¬ 
mal  management  to  constrain  the  operating  temperature  range. 
In  applications  with  high  power  demands,  such  as  automotive 
traction  batteries,  the  internal  temperature  of  the  battery  may  rise 
quickly,  due  to  joule  heating,  and  can  be  higher  than  the  surface 
temperature.  The  internal  temperature  of  the  battery  should  be 
limited  to  extend  battery  life  [14];  however,  in  practice  only  the 
surface  temperature  of  the  battery  may  be  measured.  Therefore, 
a  thermal  model  is  needed  to  predict  the  missing  information 
for  use  in  the  equivalent  circuit  electrical  model  and  for  thermal 
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management. 

Coupled  electro-thermal  models  have  been  investigated  us¬ 
ing  PDE  based  electrical  models  in  [2,4, 15],  and  equivalent  cir¬ 
cuit  based  electrical  models  in  [14, 16-18].  The  thermal  mod¬ 
els  used  in  these  studies  have  either  been  complex,  or  very  sim¬ 
ple  only  capturing  the  lumped  temperature.  Complex  thermal 
models  that  capture  the  detailed  temperature  distribution  in  a  cell 
have  been  used  [15, 19,20],  but  require  a  large  amount  of  compu¬ 
tational  resources,  making  them  unsuitable  for  control  oriented 
modeling.  A  simple  thermal  model  that  predicts  the  critical  tem¬ 
perature  of  a  cylindrical  cell  is  desired,  such  as  the  two  state  ther¬ 
mal  model  that  has  been  studied  in  [21,22].  This  model  has  the 
ability  to  capture  the  core  temperature  Tc  of  a  cylindrical  cell 
which  is  greater  than  the  surface  temperature  Ts  under  high  dis¬ 
charge  rates  [23].  Lumped  temperature  models  may  not  be  able 
to  capture  this  core  temperature.  The  two  state  thermal  model 
can  be  expanded  to  a  battery  pack  configuration  to  estimate  un¬ 
measured  temperatures  as  presented  in  [22]. 

In  this  paper  an  OCV-R-RC-RC  equivalent  circuit  electrical 
model  is  integrated  with  a  two  state  thermal  model.  In  Section  2, 
the  coupling  between  the  heat  generation  and  temperature  in  the 
integrated  electro-thermal  model  is  highlighted  by  the  tempera¬ 
ture  dependence  of  the  equivalent  circuit  parameters.  In  Section 
3,  we  first  show  how  the  electric  model  can  be  parameterized 
using  a  low  current  rate  so  that  isothermal  conditions  could  be 
assumed.  The  identified  parameters  and  their  dependencies  on 
SOC,  current  direction,  and  temperature  are  examined.  Basis 
functions  are  chosen  to  represent  the  temperature  and  SOC  de¬ 
pendence  of  the  circuit  elements.  Next  the  parameters  of  the 
thermal  model  are  identified  using  the  heat  generation  calculated 
by  the  modeled  open  circuit  voltage  for  a  high  C-rate  drive  cy¬ 
cle.  Finally  in  Section  4,  the  coupled  electro-thermal  model  is 
validated  against  the  measured  terminal  voltage  and  surface  tem¬ 
perature  data  from  a  drive  cycle  experiment. 


2  BATTERY  MODEL 

In  this  section  the  electrical  and  thermal  battery  models  are 
presented.  The  model  parameter  dependencies  are  introduced, 
and  an  electro-thermal  model  is  formed  through  a  heat  generation 
term. 

2.1  Electrical  Model 

The  battery  state  of  charge  (SOC)  is  defined  by  current  inte¬ 
gration  as, 


SOC=- 


1 

3600C,/' 


(1) 


The  nominal  capacity  of  the  cell  Cn(Ah )  is  found  by  cycling  the 
battery  cell  per  manufacturer  recommendation  [24].  The  charg¬ 
ing  profile  consists  of  a  Constant  Current  -  Constant  Voltage 
(CC-CV)  charging  cycle  that  is  terminated  when  the  current  ta¬ 
pers  below  50mA,  and  the  voltage  at  the  end  of  discharge  is  2.0  V. 


Figure  1 .  SINGLE  CELL  EQUIVALENT  CIRCUIT  MODEL. 

The  battery  electrical  dynamics  are  modeled  by  an  equivalent  cir¬ 
cuit  as  seen  in  Fig.  1.  The  double  RC  model  structure  is  a  good 
choice  for  this  battery  chemistry,  as  shown  in  [25].  The  two  RC 
pairs  represent  a  slow  and  fast  time  constant  for  the  voltage  re¬ 
covery  as  shown  by. 
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The  states  V\  and  V2  are  the  capacitor  voltages.  The  parameters 
Ri(£l),Ci(F)  correspond  to  the  first  RC  pair,  and  RiiLl)  .CMP) 
to  the  second  RC  pair.  The  states  of  the  electrical  model  are  SOC , 
Vj,  and  V2.  The  current  I  is  the  input,  and  the  model  output  is  the 
battery  terminal  voltage  Vj  defined  as, 


Vt  =  Vocv-V1-V2-IRs,  (3) 

where  Vocv  represents  the  open  circuit  voltage,  and  Rs  represents 
the  internal  resistance  of  the  cell.  The  Vocv  curve  is  assumed  to 
be  the  average  of  the  charge  and  discharge  curves  taken  at  very 
low  current  (C/20),  since  the  LiFePO 4  cell  chemistry  is  known 
to  yield  a  hysteresis  effect  as  shown  in  [25,26].  This  phenomena 
has  been  modeled  for  NiMH  and  lithium  ion  cells  [25-29],  but 
will  be  neglected  in  this  study.  The  open  circuit  voltage  Vocv 
depends  only  on  SOC;  however,  the  equivalent  circuit  parameters 
depend  on  SOC,  temperature,  and  current  direction  as  shown  in 
[10]  and  the  results  of  this  paper. 

The  cell  temperature  is  driven  by  heat  generation  Q(W )  de¬ 
fined  as, 


Q  =  I{Vocv~Vt).  (4) 

The  heat  generation  Q  in  the  battery  cell  is  defined  by  the  polar¬ 
ization  heat  from  joule  heating  and  energy  dissipated  in  the  elec¬ 
trode  over-potentials  [19].  The  effect  of  the  entropic  heat  gener¬ 
ation  is  excluded  for  simplicity,  as  it  is  relatively  small  compared 
to  the  total  heat  generation  for  an  LiFePO 4  cell  as  shown  by  [23], 
The  entropic  heat  would  contribute  less  than  1%  of  mean  Q  for 
the  drive  cycle  used  in  this  paper. 
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2.2  Thermal  Model 
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Figure  2.  SINGLE  CELL  LUMPED  THERMAL  MODEL. 


Figure  3.  ELECTRO-THERMAL  MODEL  DIAGRAM. 


The  two  state  thermal  model  is  defined  by, 

Ts  -  Tc 
CcTc  =  Q+  - - - 

Rc 

_Tf-Ts  Ts-Tc,  (5) 

S  S  D  p 

Kli  T\C 

where  Tc(°C)  and  TS(°C)  represent  the  core  and  surface  temper¬ 
ature  states  respectively.  The  temperature  used  by  the  equivalent 
circuit  model  is  the  mean  of  the  core  and  surface  temperatures 
defined  as  Tm(°C).  The  inputs  are  the  inlet  air  coolant  tempera¬ 
ture  Tf(°C)  and  the  heat  generation  Q  calculated  by  the  electrical 
model  shown  by  Eq.  4.  The  parameters  Cc(J/K)  and  CS(J/K) 
are  the  lumped  heat  capacities  of  the  core  and  surface  respec¬ 
tively,  RC(K/W)  is  the  equivalent  conduction  resistance  between 
the  core  and  surface  of  the  cell,  and  R„(K/W)  is  the  equivalent 
convection  resistance  around  the  cell.  The  convective  resistance 
R„  depends  on  the  flow  condition,  and  can  be  modeled  for  differ¬ 
ent  types  of  coolants  as  described  in  [30, 31]. 

2.3  Model  Coupling 

The  electro-thermal  model  is  formed  by  taking  the  calcu¬ 
lated  heat  generation  from  the  electrical  model  as  an  input  to 
the  thermal  model.  The  thermal  model  then  generates  the  bat¬ 
tery  surface  and  core  temperatures,  used  to  find  the  mean  battery 
temperature  for  the  parameters  of  the  electrical  model,  as  shown 
in  Fig.  3.  The  inputs  of  the  electro-thermal  model  are  the  current 
I  for  the  electrical  model,  and  the  air  inlet  temperature  Tf  for 
the  thermal  model.  The  electro-thermal  model  outputs  are  SOC, 
voltage,  and  the  battery  temperatures. 

The  coupling  results  in  a  stable  system  because  of  the  nega¬ 
tive  feedback  for  the  temperature  states.  The  negative  feedback 
coupling  between  the  thermal  and  electrical  models  can  be  ex¬ 
plained  by  the  temperature  dependence  of  battery  internal  resis¬ 
tance.  To  understand  this  coupling,  consider  a  constant  current. 


Under  this  condition  the  heat  generation  Q  will  decrease  with  in¬ 
creasing  cell  temperature,  because  the  reaction  kinetics  become 
more  favorable,  which  further  reduces  the  internal  resistance. 


3  MODEL  PARAMETER  IDENTIFICATION 

In  this  section  the  electrical  and  thermal  model  parame¬ 
terization  methods  are  described.  First  the  parameters  of  the 
equivalent  circuit  model  are  identified  from  pulse  current  dis¬ 
charge/charge  and  relaxation  experiments  at  different  SOC’s  and 
temperatures  with  the  battery  placed  inside  a  thermal  chamber. 
Then  using  the  calculated  heat  generation  of  the  cell,  the  parame¬ 
terization  of  the  thermal  model  is  presented  using  a  least-squared 
fitting  algorithm  originally  developed  in  [32], 

There  are  different  methods  of  identifying  equivalent  circuit 
model  parameters  such  as  electro  impedance  spectroscopy  (EIS) 
[11],  genetic  algorithm  (GA)  optimization  [25],  and  nonlinear 
least  squares  curve  fitting  techniques  [10].  Most  of  these  involve 
identifying  parameters  with  respect  to  SOC  as  in  [7-9,  11,33], 
in  addition  the  parameters  are  shown  to  depend  on  temperature 
and  current  direction  [10, 14,25].  The  method  selected  here  is  to 
identify  parameters  from  experimental  pulse  current  data  using 
nonlinear  least  squares  curve  fitting.  Assuming  isothermal  con¬ 
ditions  the  identification  is  performed  at  each  temperature  and 
SOC  grid  point  (5  parameters  per  pulse)  in  order  to  avoid  simul¬ 
taneous  identification  of  the  full  parameter  set  (ie.  360  parame¬ 
ters  in  this  model).  This  reduces  the  computational  burden  and 
allows  us  to  investigate  the  equivalent  circuit’s  parameter  depen¬ 
dence  on  temperature  and  SOC. 

3.1  Electrical  Model  Parameterization 

Experiments  to  parameterize  the  electrical  model  for  a 
2.3Ah  A123  26650  LiFePO 4  cell  were  conducted  using  a  Yoko- 
gawa  GS-610  Source  Measure  Unit  to  control  the  current, 
and  a  Cincinnati  Sub-Zero  ZPHS16-3.5-SCT/AC  environmental 
chamber  to  regulate  the  air  coolant  temperature.  The  tests  were 
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Discharge 


Figure  4.  PULSE  DISCHARGE  VOLTAGE  AND  CURRENT  PROFILE. 


in  the  inset  of  the  top  subplot  of  Fig.  4)  defined  as. 


Rs  = 


AK 

I 


(6) 


Where  I  is  the  current  applied  during  the  pulse  discharge/charge 
before  the  relaxation  period  (eg.  2.3A  as  shown  in  the  inset  of 
the  bottom  subplot  of  Fig.  4).  The  remaining  equivalent  circuit 
model  parameters  are  identified  by  minimizing  the  error  in  volt¬ 
age  between  the  model  and  data  during  the  relaxation  period. 


n 

J. Electrical  =  mill  iVrelaxi})  V'r.dala  (())  i  (7) 

i=l 

using  the  lsqcurvefit  function  in  MATLAB.  Each  instance  is  rep¬ 
resented  by  i,  starting  from  the  first  voltage  relaxation  datapoint 
i  =  1,  up  to  the  last  datapoint  i  =  n. 

The  voltage  recovery  during  relaxation,  Vreiax{tr),  is  derived 
by  solving  Eq.  (2),  assuming  the  capacitor  voltages  Vi,  Vi  at  the 
end  of  the  previous  rest  period  are  zero 


conducted  in  the  environmental  chamber.  The  battery  temper¬ 
ature  is  assumed  to  be  isothermal  and  Tm  equal  to  the  ambient 
temperature  in  the  chamber  due  to  the  low  C-rate  experiments. 
This  assumption  is  consistent  with  the  small  measured  rise  in 
surface  temperature  of  the  battery  cell,  less  than  0.7°C,  during 
the  pulsed  discharge. 

First  the  capacity  of  the  cell  is  measured  by  cycling  the  bat¬ 
tery  at  low  rate  (C/20).  The  Vocv  curve  is  assumed  to  be  the  av¬ 
erage  of  the  charge/discharge  curves  corresponding  to  the  same 
C/20  cycle  test  at  25°C.  The  effect  of  hysteresis  in  this  cell  chem¬ 
istry  results  in  a  voltage  gap  between  the  charge  and  discharge 
curves  as  explained  in  [25-28].  Since  hysteresis  is  not  being 
modeled  in  this  paper,  the  average  curve  is  used  for  Vocv-  It  is 
shown  in  [10],  that  there  is  a  minimal  effect  on  Vocv  with  re¬ 
spect  to  the  temperature  range  of  study  here  for  an  LiFePO^  cell. 
Therefore,  Vocv  is  modeled  with  an  SOC  dependence. 

After  Vocv  and  capacity  are  determined,  the  experiments  to 
generate  data  for  parameterization  of  the  RC  elements  are  con¬ 
ducted.  First  the  cell  sits  at  a  constant  temperature  set  point  for 
2h  to  ensure  thermal  equilibrium.  The  battery  is  then  charged 
up  to  100%  SOC  using  a  1C  CCCV  charge  protocol  at  the  3.6V 
maximum  until  a  50mA  CV  cutoff  current  is  reached.  It  is  then 
discharged  by  10%SOC  at  1C  rate,  and  relaxed  for  2h.  This  pro¬ 
cess  is  repeated  until  the  2V  minimum  is  reached.  The  pulse 
current  followed  by  a  2h  relaxation  profile  is  repeated  for  the 
charge  direction  up  to  the  3.6V  maximum.  The  pulse  discharg¬ 
ing  and  charging  is  conducted  at  different  temperatures,  resulting 
in  150C,25°C,35°C,45"C  datasets.  The  voltage  and  current  pro¬ 
file  of  one  of  the  pulse  discharge  tests  at  15°C  is  shown  in  Fig.  4. 

The  equivalent  circuit  parameter  Rs(£l)  is  found  using 
Ohm’s  law  and  the  measured  initial  voltage  jump  AVS  (shown 


Vrelax{tr)=IRl(l-exp(-t^^))(l-exp{--^—)) 

R 1C1  R\L\  ^ 

+IR2{l-exp(-  1  -  exp(- ))  +IRS . 

K2C2  A2C2 

Where  tpuise  is  duration  of  the  constant  current  pulse  prior  to  the 
relaxation  period,  and  tr  is  the  time  since  the  start  of  relaxation, 
as  shown  in  Fig.  4. 

The  inclusion  of  two  or  more  RC  pairs  in  the  equivalent  cir¬ 
cuit  model  increases  the  accuracy  of  the  cell  voltage  dynamic 
prediction  as  seen  in  [7,9, 11,25].  A  comparison  of  the  perfor¬ 
mance  for  best  fit  single  RC,  double  RC,  and  triple  RC  models  is 
shown  Fig.  5.  One  can  see  that  the  single  RC  pair  model  yields 
large  error  especially  during  the  first  500  seconds  of  relaxation, 
whereas  the  double  RC  and  triple  RC  pair  models  yield  less  error 
across  the  entire  dataset  time  period.  It  is  evident  that  the  higher 
order  RC  models  can  achieve  a  better  fit  to  the  relaxation  voltage 
data  than  that  of  the  single  RC  pair  model.  The  double  RC  pair 
model  has  sufficient  accuracy  with  the  lowest  number  of  states 
and  therefore  is  an  appropriate  choice. 

3.2  Equivalent  Circuit  Parameters 

The  equivalent  circuit  parameters  can  then  be  character¬ 
ized  as  functions  of  SOC,  and  temperature  for  the  discharge  and 
charge  direction  as  shown  in  [10].  The  calculated  internal  resis¬ 
tance  Rs  from  Eq.  (6),  is  shown  in  Fig.  6  with  respect  to  SOC  and 
temperature  for  discharge  and  charge.  The  internal  resistance  Rs 
has  a  minimal  dependence  on  SOC  over  the  range  of  10  to  90 
%,  but  depends  strongly  on  temperature  and  current  direction. 
Therefore,  the  Rs  parameter  can  be  represented  by  an  exponen¬ 
tial  function  of  the  mean  temperature  Tm,  for  the  discharging  and 
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Relaxation  after  CC  Discharge 


Figure  5.  FITTING  OF  VOLTAGE  RELAXATION  DATA. 


charging  cases,  as  shown  by. 


_jRsd,i>=  0  ( discharge ) 
RSc ,  I  <  0  ( charge ) 

TrefRs , 


Rs  = 

Ry,  =Rs0,exp( 


(9) 


T)n  TshiftRs 


) 


Where  *  =  d,c  represents  the  value  during  discharging  and 
charging  respectively.  The  characterized  Rs  functions  in  Eq.  (9) 
are  plotted  along  with  the  Rs  values  fit  from  the  relaxation  data 
using  Eq.  (7),  in  Fig.  6.  The  values  for  Eq.  (9)  are  shown  in 


Charge 


SOC 


Table  1 .  PARAMETRIC  Rs  FUNCTION  PARAMETERS. 


Rs0d 

Rs0c 

TrefRsd 

TrefRsc 

TshiftRsd 

TshiftRsd 

0.0048 

0.0055 

31.0494 

22.2477 

-15.3253 

-11.5943 

Eq.  (19)  and  Eq.  (20).  The  characterized  functions  are  plotted 
along  with  the  R i , R2  parameter  values  in  Fig.  11  and  Fig.  12. 

The  parameters  C\ ,  C2  are  represented  by  polynomial  SOC 
functions  including  temperature  dependence.  The  Ci,C2  func¬ 
tions  including  SOC  and  temperature  dependence  for  discharge 
and  charge  are  shown  in  Eq.  (21)  and  Eq.  (22).  They  are  plotted 
along  with  the  Ci,C2  parameter  values  in  Fig.  13  and  Fig.  14. 

3.3  Thermal  Model  Parameterization 

The  experiment  procedure  used  to  identify  the  thermal 
model  parameters  is  the  Urban  Assault  Cycle  (UAC),  scaled  for 
the  A123  26650  cell,  as  explained  in  [32],  This  cycle  has  been 
presented  in  [34],  for  a  13.4  ton  armored  military  vehicle.  The 
cell  is  first  charged  to  100%SOC  using  a  1C  CC-CV  protocol  un¬ 
til  the  50mA  CV  cutoff  current  is  reached.  It  is  then  discharged 
at  1C  to  about  50%SOC.  The  UAC  current  profile  is  then  applied 
to  the  cell  under  constant  coolant  flow,  with  the  measured  inlet 
temperature  Tf  as  shown  in  Fig.  7. 


Figure  6.  CALCULATED  Rs  VERSUS  PARAMETRIC  FUNCTIONS  DE¬ 
SCRIBING  THEIR  DEPENDENCE  ON  TEMPERATURE  AND  SOC. 

Tab.  1. 

The  parameters  R\.R2  are  characterized  by  including  an 
SOC  dependency  to  the  function  in  Eq.  (9)  used  for  the  param¬ 
eter  Rs.  The  corresponding  R  \ .  /C  functions  including  SOC  and 
temperature  dependence  for  discharge  and  charge  are  shown  in 


Figure  7.  UAC  CURRENT  AND  INLET  TEMPERATURE  PROFILE. 

The  experiment  is  done  by  using  a  Bitrode  FTV 1-200/50/2- 
60.  The  battery  cell  is  placed  in  a  designed  flow  chamber  as 
shown  in  Fig.  8,  where  a  Pulse  Width  Modulated  (PWM)  fan  is 
mounted  at  the  end  to  regulate  the  air  flow  rate  around  the  cell. 
This  flow  chamber  is  used  to  emulate  cooling  conditions  of  a  cell 
in  a  pack,  where  the  flowrate  is  adjustable.  Two  T-type  thermo¬ 
couples  are  used  for  temperature  measurement,  one  attached  to 
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Flow  Chamber 


Outlet 


Fan 


Inlet 


where 


z(0  = 
)  = 


r  Z(0) 

ir 

m(h) 

m(t2) 

m(t) 

¥{t2) 

m{t\) 

m(t2) 

m(t) 

\]\+¥ 

"m*)- 

(13) 


m(t)  is  the  normalization  factor  to  enhance  the  robustness  of  pa¬ 
rameter  identification  as  explained  in  [22].  For  this  purpose,  the 
parametric  model  for  the  linear  model  identification  with  initial 
battery  surface  temperature  condition  Ts  q  is  first  derived.  The 
thermal  model  in  Eq.  (5)  becomes  [22], 


Figure  8.  SINGLE  CELL  FLOW  CHAMBER. 


the  aluminum  casing  of  the  cell  to  measure  the  surface  tempera¬ 
ture  Ts,  and  the  other  near  the  battery  inside  the  flow  chamber  to 
measure  the  air  flow  temperature  Tf.  This  thermal  identification 
experiment  setup  is  also  presented  in  [32]. 

The  non-recursive  least  squares  thermal  model  identification 
method  described  in  [22,  32]  is  implemented  here  by  using  the 
heat  generation  from  Eq.  (4)  as  the  input  for  the  thermal  model. 
Where  Vj  is  the  measured  voltage,  Vocv  is  the  modeled  open 
circuit  voltage,  and  I  is  the  measured  current.  The  objective  is 
to  minimize  the  sum  of  the  squared  errors  between  the  the  mod¬ 
eled  surface  temperature  Ts,  and  the  measured  surface  tempera¬ 
ture  Ts  i]ala  as  shown  by  the  cost  function. 


t  O - L  "'I,  J  L - U 

{  CC  T  CS  |  1  rp  T  \ 

—  \rr  n  +  7^~R~^STs  ~  ^>0' 


(14) 


after  a  Laplace  transformation  and  substitution  of  the  unmeasur¬ 
able  Tc  by  the  measurable  7)  ,  Ts.  In  order  for  this  algorithm  to  be 
implemented,  a  proper  parametric  model  must  be  used.  A  filter 
is  designed  and  applied  to  the  parametric  model  in  Eq.  (11), 


(15) 


where  the  observation  z  and  the  independent  regressors  (f>  are 
measured.  The  parameter  vector  0  is  defined  as  [22], 


n 

^Thermal  =  min  ^  (7j(0  ^s.dala  (f) )  ■  (10) 

i=l 


Where  each  instance  is  represented  by  i,  starting  from  the  first 
surface  temperature  datapoint  i  =  1,  up  to  the  last  datapoint  i  =  n. 
A  parametric  model  in  the  form  of  [35], 


z  =  s2Ts  -  sTsfi 
ty=[QTf-TssTs-Tsfl]T 

e  =  [«  P  y]  t 


where  the  parameters  a,  p,y  are  [22], 


(16) 


z  =  er4, 


(ii) 


is  used  for  the  thermal  model  parameter  least  squares  identifica¬ 
tion  [22].  Where  the  observation  z  and  the  independent  regres¬ 
sors  (f)  should  be  measured.  The  parameters  in  0  are  calculated 
by  the  non-recursive  least  squares  after  the  experimental  data  is 
taken  over  a  period  of  time  ...,f  by  [35], 


1 


a  = 


P  = 


CcCsRc 

1 

CcCsRcRu 


Y  = 


CcCs 

CcCsRc 


1 


CSR 


U 


) 


(17) 


0(o  =  (*i’(f)*(or1«‘(oz(o, 


Assuming  Cc,  the  parameters  a.  p.y  are  identified  to  find  CS,RU, 
(12)  and  Rc.  This  is  achieved  based  on  the  measured  observation  and 


6 


UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release 


Table  2.  THERMAL  MODEL  PARAMETERS. 


CS{J/K) 

Cc{J/K) 

RC(K/W ) 

RU(K/W) 

8.3 

71.95 

1.1139 

3.1421 

regressor  vectors  shown  as  [22], 

(Cc  +  Cs)CsfiR„2  +  CsjRu  +  1=0 

1 

Rc  "  pCsCcRu 

Cs  =  — - — 

CcaRc 


(RMSE)  in  predicted  surface  temperature  is  0.33 °C  and  voltage 
is  18.4mV.  The  voltage  RMSE  is  comparable  with  published  re¬ 
sults  in  [29],  using  a  similar  type  of  drive  cycle  profile  for  this 
type  of  cell.  The  predicted  core  temperature  Tc  is  also  shown  in 
Fig.  10,  which  is  1.77°C  higher  than  the  predicted  surface  tem¬ 
perature  Ts  under  this  cycle.  A  higher  7]  prediction  is  presented 
in  [32],  using  a  different  heat  generation  under  the  same  exper¬ 
imental  conditions.  The  heat  generation  for  our  case  is  smaller 
than  [32],  causing  slightly  different  identified  parameters  and  a 
lower  Tc  prediction.  Further  investigation  is  required  to  deter¬ 
mine  if  the  calculated  heat  generation  Q  and  core  temperature 
T(  prediction  are  correct.  Including  a  hysteresis  model  in  the 
electro-thermal  model  will  also  need  to  be  investigated  to  deter¬ 
mine  if  better  results  can  be  achieved. 


The  parameter  Cc  is  assumed  to  be  close  to  that  in  [23,  32] 
and  the  resulting  identified  parameters  CSlRc  and  R„  from  the 
thermal  identification  scheme  are  shown  in  Tab.  2.  The  parame¬ 
ters  Cc,Cs,Rc  should  not  change  significantly  within  the  lifetime 
of  the  battery  cell  due  to  their  physical  properties.  The  parameter 
Ru  can  change  with  respect  to  the  flow  around  the  cell  as  previ¬ 
ously  mentioned.  In  this  case  it  is  identified  as  a  constant  for  a 
steady  flow  condition. 

4  MODEL  VALIDATION  AND  RESULTS 

The  electro-thermal  model  is  implemented  in  Simulink  to 
validate  its  performance  under  the  UAC  experiment.  The  SOC, 
temperature,  and  current  direction  dependencies  of  the  equiva¬ 
lent  circuit  model  parameters  are  included  using  lookup  tables. 
The  current  I  and  air  inlet  temperature  Tf  inputs  are  shown  in 
Fig.  7.  The  voltage  and  temperature  responses  of  the  electro¬ 
thermal  model  are  compared  to  the  experimental  measurements. 
SOC  is  shown  in  Fig.  9  for  this  experiment. 


Figure  9.  UAC  SIMULATION  SOC  RESULTS. 


The  SOC  varies  between  50%  and  43%  under  these  conditions. 

The  measured  surface  temperature  Tsdata  and  terminal  volt¬ 
age  Vr^data,  are  compared  to  the  predicted  surface  temperature  Ts 
and  voltage  Vj  as  shown  in  Fig.  10.  The  root  mean  square  error 


_ T _ T  _ T  .  , 

c  s  s,data 


_ I _ | _ | _ | _ | _ i  OU  QU.U  Oil 

0  5  10  15  20  25  30  35  40 


Time  (min) 

Figure  10.  UAC  TEMPERATURE  AND  VOLTAGE  RESULTS. 


5  CONCLUSION  AND  FUTURE  WORK 

In  this  study  an  equivalent  circuit  electrical  model  along 
with  a  two  state  thermal  model  for  an  A 123  26650  LiFePO 4 
cell  were  parameterized.  The  models  were  integrated  into  an 
electro-thermal  model  in  MATLAB/Simulink  through  a  coupling 
heat  generation  and  temperature  feedback.  The  resulting  electro¬ 
thermal  model  matches  experimental  measurements  with  mini¬ 
mal  error.  This  shows  that  the  parameterization  schemes  used 
are  adequate  for  battery  modeling. 

Future  work  will  involve  modeling  of  hysteresis  as  in  [25, 
27,28],  which  will  then  cause  the  heat  generation  to  change  due 
to  the  new  Vocv  term.  Measurement  of  the  core  temperature  Tc 
is  also  planned  to  validate  the  electro-thermal  model. 
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Appendix:  Calculated  Parameters  and  Functions 
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Ru  =(Riot+Rnt(SOC)+Rl2t(SOC)2) 


exp{ 


*  refR  i  * 


Tm  TshiftR\* 


) 


(19) 


Discharge 


SOC 


Charge 


• 

15°C 

— 

15°C 

Fit 

25°C 

— 

25°C 

Fit 

— 9— 

35°C 

— 

35°C 

Fit 

X 

45°C 

— 

45°C 

Fit 

Figure  11.  CALCULATED  R\  VERSUS  PARAMETRIC  FUNCTIONS 
DESCRIBING  THEIR  DEPENDENCE  ON  TEMPERATURE  AND  SOC. 


Table  3.  PARAMETRIC  R\  FUNCTION  PARAMETERS. 


Riod 

R  i  oc 

R\\d 

Rnc 

R\id 

7.1 135e-4 

0.0016 

-4.3865e-4 

-0.0032 

2.3788e-4 

R\lc 

TrefR\d 

TrefR\c 

Ts  hi ft  R\d 

TshiftR\c 

0.0045 

347.4707 

159.2819 

-79.5816 

-41.4548 

_  (R2d,  I  >=  0  ( discharge ) \ 

2  \  R2c  ,  /  <  0  ( charge )  J 

R2,  =  (Rio,  +R2h(SOC)+R22t(SOC)2)eXp(7^^)  (20) 


Figure  12.  CALCULATED  R2  VERSUS  PARAMETRIC  FUNCTIONS 
DESCRIBING  THEIR  DEPENDENCE  ON  TEMPERATURE  AND  SOC. 
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Table  4.  PARAMETRIC  R2  FUNCTION  PARAMETERS. 


R20d 

R20c 

R2ld 

R2\c 

f  C2rf ,  I  >=  0  ( discharge )  1 

2  [C2c,  I  <  0  (charge)  J 

0.0288 

0.0113 

-0.073 

-0.027 

R22a 

R22c 

T'refR'id. 

TrefR2C 

C2*  =C2o* +C2i*(50C) +C22*(50C)“ 

(22) 

0.0605 

0.0339 

16.6712 

17.0224 

+  (C23*  +  C24*  (SOC)  +C25*(SOC)~)Tm 

r_[  C\d ,  I  >=  0  ( discharge )  1 
1  \Cic,  /<  0  ( charge )  J 

CU  =Cio*  +  Cn*(SOC)  +Ci2*(SOC)~  (21) 

+  (C13*  +  C14*  (SOC)  -\-Ci5*(SOC)~)Tm 


1  — 15°C 

- 15°C 

— ^ — 25°C 

- 25°C 

— B—  35°C 

- 35°C 

Fit 

X  45°C 

|  - 45°C 

Fit 

soc  soc 

Figure  13.  CALCULATED  Cj  VERSUS  PARAMETRIC  FUNCTIONS 
DESCRIBING  THEIR  DEPENDENCE  ON  TEMPERATURE  AND  SOC. 


Table  5.  PARAMETRIC  C\  FUNCTION  PARAMETERS. 


Ci0d 

ClOc 

Cud 

Cm, 

335.4518 

523.215 

3.1712e+3 

6.4171e+3 

C\2d 

C12c 

Cl3d 

Cl3, 

-1.3214e+3 

-7.5555e+3 

53.2138 

50.7107 

Cud 

Cuc 

Cl5rf 

Cl  5, 

-65.4786 

-131.2298 

44.3761 

162.4688 

Figure  14.  CALCULATED  C2  VERSUS  PARAMETRIC  FUNCTIONS 
DESCRIBING  THEIR  DEPENDENCE  ON  TEMPERATURE  AND  SOC. 


Table  6.  PARAMETRIC  C2  FUNCTION  PARAMETERS. 


C20d 

C20c 

C2irf 

C2i, 

3.1887e+4 

6.2449e+4 

-1.1593e+5 

-1.055e+5 

C22d 

C22, 

c23rf 

C23, 

1.0493e+5 

4.4432e+4 

60.3114 

198.9753 

C24d 

C24, 

C25  d 

C25, 

1.0175e+4 

7.5621e+3 

-9.5924e+3 

-6.9365e+3 
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